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We present a theoretical study of gap opening in the zeroth Landau level in gapped graphene as a 
result of pseudo-Zeeman interaction. The applied magnetic field couples with the valley pseudospin 
degree of freedom of the charge carriers leading to the pseudo-Zeeman interaction. To investigate 
its role in transport at the Charge Neutrality Point (CNP), we study the integer quantum Hall 
effect (QHE) in gapped graphene in an angular magnetic field in the presence of pseudo-Zeeman 
interaction. Analytical expressions are derived for the Hall conductivity using Kubo-Greenwood 
formula. We also determine the longitudinal conductivity for elastic impurity scattering in the first 
Born approximation. We show that pseudo-Zeeman splitting leads to a minimum in the collisional 
conductivity at high magnetic fields and a zero plateau in the Hall conductivity. Evidence for 
activated transport at CNP is found from the temperature dependence of the collisional conductivity. 
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I. INTRODUCTION 

In recent years, the experimental realization of a stable 
single layer of carbon atoms[l, 2] has stimulated much in- 
terest in the studies of its unusual properties [3, 4]. This 
material known as graphene is a tightly packed honey- 
comb lattice of carbon atoms. Graphene monolayer is a 
gapless semiconductor with conical touching of electron 
and hole bands. The charge carriers in this system obey a 
linear dispersion relation near the Dirac point, which en- 
dows it with unique electronic properties. This difference 
in the nature of the quasiparticlcs in graphene from con- 
ventional two-dimensional (2D) electronic systems has 
given rise to a host of new and unusual phenomena. Be- 
sides the fundamental interest in understanding the elec- 
tronic properties of graphene there are also serious efforts 
to build nanoelectronic devices from graphene [5-7]. 

The quantum Hall measurements in graphene were one 
of the key tools providing evidence that the quasiparti- 
cles in graphene are chiral, massless fermions known as 
Dirac fermions [8-13]. A key difference in the integral 
quantum Hall effect in graphene compared to the stan- 
dard effect in conventional two-dimensional electron gas 
(2DEG) systems is the occurrence of a zeroth Landau 
level (LL) state. There have been several experimental as 
well as theoretical studies of transport at the CNP which 
required analyzing the role of the zeroth LL state[14-19]. 
The nature of the splitting of the electronic states at the 
CNP at present remains unclear, whether this splitting is 
due to Zeeman interactions, electron-electron or electron- 
phonon interactions, asymmetric gap, gapless edge states 
or due to valley splitting is still an open question. The re- 
sults obtained in [16] were explained on the basis of field 
dependent splitting of the zeroth Landau level. Break- 
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down of the QHE in graphene leading to two insulating 
regimes were analyzed in [17, 18]. On the theoretical 
side, [15] addressed the role of the zeroth LL in the QHE 
in graphene using Laughlin's gauge argument. The role 
of disorder in addition to the Landau level structure was 
discussed in detail in [19]. In [14], the splitting of the 
zeroth Landau level was observed and it was attributed 
to lifting of the sublattice and spin symmetry. 

We consider gapped graphene in the presence of an ex- 
ternal magnetic field. Due to the breaking of sublattice 
symmetry a gap opens in the energy spectrum at the 
Dirac points. One of the mechanisms that can lead to 
sublattice symmetry breaking is through a local asym- 
metric chemical or electrical environment provided by 
a substrate, such as epitaxial graphene on SiC or BN 
substrate. Band gaps of various magnitudes can be in- 
duced depending on the substrate [20, 21]. This symme- 
try breaking can have important consequences for trans- 
port at CNP as it can lead to pseudo-Zeeman splitting of 
the zeroth level when contributions of both the valleys are 
considered. It has been shown in [23, 30-33, 40, 43, 45] 
that if Berry phase effects are taken into account, for 
crystals with broken spatial inversion symmetry, the elec- 
trons acquire an orbital magnetic moment as a result 
of the self rotation of the Bloch electron wave packet. 
This applies to Dirac fermions in graphene with stag- 
gered sublattice potential which breaks inversion sym- 
metry where the orbital magnetic moment is associated 
with the valley index and can lead to the valley QHE. In 
addition, various mechanism for valley filtering and valley 
polarization have been discussed in relation to electronic 
devices[33, 40-43]. 

Our focus is on electron transport at CNP in gapped 
graphene in a tilted magnetic field when the sublattice 
symmetry of graphene is broken resulting in an energy 
gap at the Dirac point. This requires that we con- 
sider valley QHE which has been discussed earlier in 
[22, 33, 40] and more recently in [23]. In [22], the authors 
relate the heights of the plateaux in the Hall conductivity 
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and the peaks in the diagonal conductivity to the size of 
the bandgap and the amount of disorder in the system. A 
semiclassical presentation of the valley QHE in graphcne 
is given in [23] . Here, we present a full quantum mechan- 
ical transport theory for valley QHE and analyze the ef- 
fects of the external tilted magnetic field on the transport 
at CNP at finite temperature in the presence of screened 
charged impurities. From the very beginning of our cal- 
culation, we explicitly introduce the Zeeman coupling of 
the external magnetic field with the valley pseudospin of 
the Dirac fermions and diagonalize the graphene Hamil- 
tonian in its presence. Moreover, we consider an external 
angular magnetic field [14, 15] applied to the system in or- 
der to highlight the role of the out-of-plane to in-plane 
magnetic field in the splitting of the zeroth LL. Further- 
more, the analysis of magnetotransport in this work is 
performed in the presence of elastic scattering due to 
charged impurities which are known to be the dominant 
scattering mechanism[24-29] in graphene on a substrate. 
In addition, this is a finite temperature study where the 
role of temperature in magnetotransport at CNP can be 
investigated. The electrical transport coefficients have 
been obtained from the standard Kubo formula in the 
self-consistent Born approximation^, 9, 30]. 

In section II, we present the formulation of the problem 
and numerical discussion of density of states for different 
magnetic field strengths with varying tilt angle from out- 
of-plane to in-plane. Section III contains the derivation of 
the Hall conductivity as a function of the tilted magnetic 
field including discussion of numerical results where as 
in section IV we evaluate the longitudinal conductivity 
as a function of tilt angle with discussion of results. In 
section V, summary of the work is given followed by an 
appendix in section VI. 



II. FORMULATION 

We consider Dirac fermions in graphene which is in the 
x — j/-plane in the presence of a tilted magnetic field and 
pseudo-Zeeman interaction. The magnetic field {{B x , 0, 



B z ) = (Bs'md, 0, Bcos9)} is applied at an angle 9 with 
the z-direction which is perpendicular to the graphene 
plane. The effective Hamiltonian for Dirac fermions in 
gapped graphene in a magnetic field[31-33] can be ex- 
pressed as (the speed of light c = 1 in the minimal sub- 
stitution that follows) 

H T * = V F [a x (p x + eA x )T z + <7 y (p y + eA y )] + A z a z . (1) 

Here t z = ±1 for valleys K and K', (A x ,A y ) are the 
components of the vector potential, Vp characterizes the 
Fermi velocity of Dirac fermions. We identify A z = 
—/-i* B B z as the pseudo-Zeeman term with B z = Bcos9, 
& = {^x^y^z} are the Pauli matrices, the effective 
Bohr magneton is (i* B = with the effective mass 

m* = 2Ah 2 /3a 2 t 2 . The Bohr magneton and the effective 
mass are expressed in terms of the gap energy A = 0.28 
eV (for graphene on SiC), the nearest neighbour hopping 
energy t = 2.82 eV and the lattice constant a = 0.246 nm 
with the result that the effective Bohr magneton (n%) can 
be 30 times larger than the free electron Bohr magneton 
(hb) [23, 33, 40]. This also allows us to ignore the real 
spin Zeeman term. The above Hamiltonian is the same as 
the one obtained in [45, 46] in the absence of valley-orbit 
coupling. 

The Hamiltonian H Tz for the two valleys (K, K') can 
be written as 



ff T * = V F 



PxTz 



A Z /V F 

+ ip y + ieAy 



PxT z - ipy ~ ieA y 
-A Z /V F 



(2) 



where in the diagonal terms A z = ±fi* B B z represents the 
potential asymmetry between A and B lattice sites, which 
opens an energy gap at the CNP. We have employed the 
Landau gauge and expressed the vector potential as ~A = 
(0, B z x — B x z, 0). The last term in the Hamiltonian 
given in Eq.(l) is regarded as the pseudospin Zeeman 
term (A 2 = fi* B B z ){23, 30-33, 43], where the valleys K 
and K' serve as pseudospin up (+1) and pseudospin ( — 1), 
respectively. Equation (2) is expressed as 



J 



H T '~ = V F 



A Z /V F p x r z 
p x r z + ip y + ie(B z x - B x z) 



iPy 



ie(B z x - B x z) 
-A Z /V F 



r 



(3) 



To obtain the energy eigen solutions of the above equa- 
tion, one can use the eigenvalue equation for a given 
spinor 



*(r) = 



Mr) 



(4) 



H 



Mr) 



E 



Mr) 



(5) 



K M r )J \Mr)y 
It yields the following equations 

A„0i(r)-iVj? (ip x r z +p y + e(B z x - B x zj) 2 (r) = S0i(r) 

(6) 
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and for n = 0, the above equation at CNP, is written as 



iV F (~ip x T z +p y + e(B z x - B x z)) 0i(r)-A^ 2 (r) = E<j) 2 (r). 

(7) 
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We try the ansatz 
Vi(r)\ 



exp[ik y y] 



V2(r) 



(8) 



where L y is the dimensions of the graphene monolayer 
in the y-dircction. From Eqs.(6,7 & 8) we obtain the 
eigenvalues of the nth LL as 



El* 



-t z A z , 



n = 



(9) 



Kn = sJ&wfr \n\ + (r 2 A 2 ) 2 , n + 



with s = ± for electrons and holes. The corresponding 
cigenfunction is 



Jk y y 



(x+xq) 



sci^|„|_il^ — r 



(10) 



The n = Landau level (the zeroth LL) requires separate 
treatment. It lies just at the top of the valence band and 
its amplitude is only at B sites with energy E$ = —A z , 
ci = and c 2 = 1. In the above equation — is for holes 

and + for electrons[31, 32], xo = l 2 k„, I 
c\ = sin(a„/2) with sin(a„) = - 
cos(a„/2) with cos(a„) 



h_ 

Y eBcost 
shuj D ^j2\n\ 

The coef- 



y/h^ D 2\n\+(T z A z y 

ficients C\ and c 2 are normalized such that \c\\ + | C2 1 2 = 
1. Vn[^] = ^k^H n (^)eM-^n 

H n (x) are the Hcrmitc polynomials, uj£> = Vf-J sM^sl. 

is the cyclotron frequency of Dirac Fermions and n is 
an integer. Similarly, for the K' valley (t z = —1), the 
Hamiltonian yields the same eigenvalues as given in Eq. 
(9) with eigenfunctions 



_ e**y ( s 'cW\ n \[^±pl] 



C2<P\n\-ll— 



(11) 



where s' = ±for electrons and holes, c[ = sin(a^/2) 



with sin(o4) 
cos(a' n ) 



y/h^l2\n\ + (T z A z y 



■■, c' 2 = cos(a^/2) with 
The energy of the n = 



^/fi 2 ^2|n| + (r z A i ) 2 

LL for the K' valley is E r ^~~ x = A z . In this case, the 
Landau level n = lies just at the bottom of the conduc- 
tion band and its amplitude is only at the A sites. The 
Density Of States (DOS) is defined as 
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(12) 



D CNP (e) = —J2He - E T >) 



where Eq z = =fA z for K and K' valleys respectively, 
and A is the area of the sample. Assuming a Gaussian 
broadening of width T, the DOS at CNP is expressed as 



D 



CNP 



2ttZ : 



1 



rV27r 



exp 



2T 2 



. (13) 



where T is the Gaussian distribution broadening width 
of zero shift. Similarly, one can evaluate the DOS for the 
other valley (K'). 
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FIG. 1: Density of States (DoS), dimensionless as a function 
of the Fermi energy for different values of the magnetic field 
strength. The magnetic field strength varied from 5 Tesla 
(solid line), 10 Tesla (dot-dashed line), 20 Tesla (dotted line) 
to 30 Tesla (dashed line) for fixed values of angle (zero degree) 
and temperature (T = 0). 



The above expression for the density of states at CNP 
is plotted in Fig. (1) as a function of the the Fermi energy 
(gate voltage) as the magnetic field strength is varied to 
see the splitting of the zeroth Landau level. The applied 
magnetic field is perpendicular to the graphene plane, 
the tilt angle 9 — 0. The magnetic field is varied from 
5 Tesla (solid line), 10 Tesla (dot-dash line), 20 Tesla 
(dotted line), all the way to 30 Tesla (dashed line). In 
Fig. (1), we see that as we increase the magnetic field 
strength, the zeroth LL splits further apart. The gap in 
the density of states is well resolved for high magnetic 
fields. 

In Fig. (2), we plot the density of states at CNP as 
a function of the Fermi energy for different values of the 
tilt angle of the magnetic field (from our-of-plane to in- 
plane). The magnetic field is fixed at 30 Tesla. The tilt 
angles are chosen to be 9 = (solid line), 40 (dot-dashed 
line), 60 (dotted), and 80 (dashed) degrees in Fig. (2). 
The splitting of the zeroth Landau level washes out as 
we increase the tilt angle of the magnetic field which is 
consistent with the pseudo-Zeeman interpretation of this 
splitting. At 9 — 80 degrees, when the magnetic field is 
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FIG. 2: Density of States (DoS), dimensionless as a func- 
tion of the Fermienergy for fixed value of the magnetic field 
strength at 30 Tesla. The tiltangle is varied from zero (solid 
line), 40 (dot-dashed line), 60 (dotted line) to 80 degree 
(dashed line) at zero temperature (T = 0). 



almost completely aligned with the graphene plane, the 
gap has closed and a single peak with no splitting occurs 
at CNP. In Figs.(l and 2) the following parameters were 
employed: g = 60[31, 32]. The broadening of the Landau 
levels generally depends on the magnetic field strength, 
the temperature, Landau level index and the scattering 
parameters. This requires a self consistent calculation 
usually performed numerically. In order to carry out a 
tractable analytical calculation we have a chosen a con- 
stant level width Y = 10 meV[25, 34-36]. 

In the presence of a magnetic field, there are two con- 
tributions to magnetoconductivity[37, 38]: the collisional 
(hopping) contribution and the Hall contribution. The 
former is the localized state contribution which carries 
the effects of SdH oscillations. The Hall contribution is 
the non diagonal contribution. In order to calculate the 
electrical conductivity in the presence of pseudo-Zeeman 
interactions and a tilted magnetic field we will follow 
the formulation of [39], which is derived from the gen- 
eral Liouville cquation[37, 38] and includes dissipative 
effects. In the linear response regime, the conductivity 
tensor is a sum of a diagonal and a nondiagonal part 



/ / V 



M + a ni(u), fi>,v = x,y. In general, 



the diagonal conductivity a^(uj) = a d ^J f (u) + ct™ 1 ^), 
accounts for both diffusive and collisional contributions 
whereas the Hall contribution is obtained from the non- 
diagonal conductivity a^(uj). Here, a x l J f = cr** ' = 
(here the diffusion contribution is zero because the di- 
agonal elements of the velocity operators vanish) and 
a x ° x = cr™ 1 . This formulation has been employed success- 
fully in electronic transport in 2DEG systems [37, 38] and 
more recently in graphene [39], and references therein. 



III. HALL CONDUCTIVITY 

The Hall conductivity a yx is obtained from the nondi- 
agonal elements of the conductivity tensor, given by [3 7- 
39] 



2ihe 2 



yx ~ n 

X 5-^ 



£/«(i-/e)<£K 10 <£'K 10 



(14) 



Since /fa-JVXe^-^) = /e(W?) and fi -> S = 
L x L y , we obtain 

<Jyx = —J— 2^ (h - J(') — ' I 15 ) 



•So 



{E 6 - E e 



Since the x and y components of the velocity operator are 
v x = |£- and v y = §S- where H T ' = V F [cr x p x T z +(7y(py + 
eA y )] + A z a z . Therefore, v x = V F a x r z and v v = V F a y . 
Hence 



V x |0 = CiC 2 VF{8 n -l, n i + S n ,n'-l) 



(16) 



and 



(£1 v y 10 = -«CiC 2 VF(£n'-l,n ~ ^n-l.n')- (17) 



Since |£) = \n,s,T z ,k y ), there will be one summation 
over k y which, with periodic boundary conditions for k y , 
will give 



dk - — 



(18) 



-L X /2P 



Substituting the values of the matrix elements of velocity 
in Eq. (15) yields 



_ 2 x 2( Cl c 2 ) 2 he 2 V F 2 
V ~ 2^P 

(/{ _ h') [$n,n'-l — S n -i t . 



(19) 



Here factor of 2 is due to spin degeneracy. Since 
W s % = s^Jh 2 uj 2 D 2 \n\ + (t z A z ) 2 we obtain 



(£ e -%r = [s^2K>w 2 D \n\ + {r z A z f 

-/^MTlW] 2 . (20) 

Substituting Eq. (20) in Eq. (19) we obtain the Hall 
conductivity 



The above equation can be further simplified and the final 
result for the angular Hall conductivity (see Appendix for 
details) is 



2(cic 2 ) 2 e 




X (f+,n ~ I- 



where we have introduced the sum over j = ±1 for 
a concise final expression. The effect of the tilted 
magnetic field can be seen in the distribution func- 

tion through oj d = V F J^f^ as /£ a = f(E^ s ) = 



r r s^/n^l2\n\+(jA z y-E f 
l eX P^ k B T 

ity at CNP is written as 
2(cic 2 



l)]" 1 . The Hall conductiv- 




n=0,j 



(/i >0 -/i,i+/i,o-£,i) 



where f> = f(E> ) = [exp( ^^ + l)]" 1 . In the 
limit when the tilt angle (6 = 0) and the Pseudo- 
Zeeman term vanishes (A z = p* B B z = 0), and we con- 
sider transport contribution from a single valley only, the 
results obtained are consistent with previous works in 
the literature[l-3, 8-10, 37]. Elements of the resistiv- 
ity tensor p flv (p,iy=x,y) can be determined from those 
of the conductivity tensor a^, obtained above, using 
the expressions: p xx = a yy /S, p yy = a xx /S and 
Pxy = -cryx /S where S = cr xx a yy - cr xy a yx with S w 
a% = n 2 e e 2 /B 2 cos 2 9. 

The Hall conductivity, as a function of the Fermi en- 
ergy, for different values of the temperature is shown in 
Fig. (3). In this figure, we have shown results for dif- 
ferent values of temperature 5 (solid line), 25 (dotted 
line), 75 (dashed) and 125 K (dot-dashed). The step 
around CNP is washed out completely at 125 K. The 
magnetic field is 5 Tesla, g = 60. This shows that in 
order to observe pscudo-Zeeman splitting of the zcroth 
LL one needs to be in the regime of high magnetic fields 
such that the thermal energy does not wash out the split- 
ting of the zeroth LL. Here, we have not considered the 
electron-phonon interaction which can affect the results 
at high temperatures. That is intended for future work. 
Furthermore, as the tilt angle of the magnetic field is 
increased, the perpendicular component of the magnetic 
field becomes smaller and the plateaus in the Hall con- 
ductivity disappear. We find that at a temperature of 
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FIG. 3: Vanishing of the pleatues in the Hall conductivity as 
a function ofthe Fermi energy with change in temperature. 
Temperature is varied from 5 (solid line), 25 (dotted line), 
75 (dashed line) to 125 K (dot-dashed line)for fixed values of 
magnetic field (5 Tesla) and the Zeeman factor (g=60). 



5K and magnetic field 5 Tesla, when 8 — 80 degrees, the 
steps disappear completely which is consistent with the 
discussion of the density of states in section II. 



IV. COLLISIONAL CONDUCTIVITY 

To obtain collisional contribution to conductivity, we 
assume that electrons are elastically scattered by ran- 
domly distributed charged impurities as it has been 
shown that charged impurities play a key role in the 
transport properties of graphene near the Dirac point. 
This type of scattering is dominant at low temperature. 
The collisional conductivity when spin degeneracy is con- 
sidered is given by [37-39] 



^^/(S)(l-/(i?'))^(^i?')(^-^) 2 



(23) 



where f{E) = [exp( 



e ~ Ef -J- 1)]" 1 is the Fermi Dirac dis- 
f(E') for elastic scattcr- 



k B T 

tribution function with f{E) 
ing, k,B is the Boltzmann constant and p is the chemical 
potential. is the transmission rate between the one- 
electron states |£) and So the volume of the system, 
and e the electron charge. Conduction occurs by transi- 
tions through spatially separated states from to xp, 
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where x^ = (£| x |£) is the mean value of the x com- 
ponent of the position operator when the electron is in 
state |£) . This is the well known hopping type formula for 
transport in the presence of a constant external magnetic 
field. Collisional conductivity arises as a result of migra- 
tion of the cyclotron orbit due to scattering by charge 
impurities. The scattering rate W^> is given by 



W (( .(E,E') 



— £ 

9 



S H 



\U q V FAu) 



5{E-E')5 ky . K+qy . 

(24) 



The Fourier transform of the screened impurity potential 
is 



U q = 



k , 



where Uo = e 2 /47reoe; ko is the screening wave vector, 
e is the static dielectric constant of the material and eo 
is the dielectric permittivity. F , (u) are the form fac- 
tors, (f| e l "- r \C) with u = l 2 {q 2 x + q 2 y )/2 = ^~ with 

2 / 9 . 9\ »x . ,i 

9± 



- (q^ + q 2 ). Ni is the impurity density and the wave- 



function is |£) = \n,s,T z ,k y ). In the situation studied 
here the diffusion contribution is zero because the diago- 
nal elements of the velocity operators vanish. Noting that 
a xx = a yy an d f° r screened impurity scattering such that 
ko » q, we can ignore the q dependence in Eq. (25). 

Here (£| x |£) = Xo is the expectation value of the po- 
sition with (x£ — x^i) 2 = (l 2 q y ) 2 and q y = q± sm(. Since 
the wave function oscillates around — xq, we have 



E, L y 
2ir 



dk h 



So 
2irP 



(26) 



-L x /2i 2 

and using cylindrical coordinates, 



E 



So 



2jt 



<i£ / du. 



(27) 



The following matrix element between the two states can 
be evaluated to yield 



(u) + c\F n 



(28) 



e u \c\L n (u) + c\L n _x (u)] 2 ; n = 



with 



Finally, evaluating the above integral, we obtain the fol- 
lowing result 



col 



where 

oo 

/ due~ 
o 



e 2 2AW 2 /? 
h PkyHu/D 

xf(EZ)(l-f(EZ)), 

we have 



^2 [c\(2n + l) + ^{2n-l) + c\c 2 2 (2n)\ 



(31) 



[c\L n (u) + cfL n . 



used 



the 

cf(2n 



relation 
+ 1) + 



!(2n 



l) + c 2 c 2 ' 



action: c% 



.^(2n). In the limit of zero Zeeman inter- 
l, 2 — t~i — j and the integral will yield 6n, which 
is consistent with previous theoretical results [8-10, 39]. 
In the above expression, the tilt angle of the magnetic 
field and the Zeeman interaction contribution appears in 

and the distribution 
I)]" 1 



(25) the magnetic length I 



■Bi 



[exp( 



k B T 



function f{E%*) 
respectively. 

At the CNP, where the contribution of the n = 
Landau level is crucial, the collisional conductivity is ex- 
pressed as 



(at CNP) cc 5>/(£^)(l - f(E?)). 



(32) 
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FIG. 4: Gap opening in the collisional conductivity as a func- 
tion ofthe Fermi energy at low temperature. Temperature is 
varied from5 K (solid line), 25 K (dotted line), 75 K (dashed 
line), to 125 K(dot-dashed line) for fixed values of the mag- 
netic field (5 Tesla)and tilt angle (0 degree). 



\F„,Au)\ 2 = — e-V 



Ll7 n (u) ;n'<n, (29) 



Inserting Eq. (24, 25, 26, 27 & 28) in Eq.(23) the colli- 
sional conductivity can be written as 



2(3Ni 



h l 2 huj 

x/(K:J(i 



E 



u- 



D 



due 



[c\L n (it) + c 2 L„_i («)] 



The collisional conductivity at CNP given by Eq. (32) 
is shown graphically in Fig (4) as a function of the Fermi 
energy for a fixed magnetic field at 9 = tilt angle as the 
temperature (T) is varied: T= 5K (solid line), 25K (dot- 
ted line),75K (dashed line) and 125K (dot-dashed line) in 
2Fig. (4). The parameters [16-1 9, 24] used in all of our fig- 



ures are : Ni = 3 x 10 



15 



S k = 10- 



V F = 10 b 



f(K%)) 



m/s, magnetic field is 5 Tesla and U = e /47re e. We 
(30)bakc e = 4 (graphene on a Si02 substrate) and eo the 
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dielectric permittivity of free space, with kp = (Trn-e) 1 ^ 2 
being the Fermi wave number. In Fig. (4), as the temper- 
ature is decreased, the collisional conductivity exhibits a 
gap around CNP due to splitting of the zeroth Landau 
level. This splitting at CNP is due to the pseudo- Zee- 
man interaction as discussed earlier in the context of the 
density of states at CNP. With an increase in tempera- 
ture, at 195 K, the splitting is completely washed out and 
there is only a single peak at CNP. Furthermore, Eq.(32) 
in the limit of low temperatures or high magnetic fields, 
yields the temperature dependence of the collisional con- 
ductivity as (T™ 1 (at CNP)oc f3e~@ Az which represents 
an activated type of behavior of the conductivity at the 
CNP. 

In the limit of zero temperature (T = 0), equation (32) 
can be expressed as 

a x°x ( at CNP) cc (e — Eq z ) , (33) 

which can then be written in terms of the Gaussian den- 
sity of states as derived and discussed in section II. The 
effect of the angular magnetic field on the collisional con- 
ductivity follows the discussion presented in section II for 
the density of states. Further, the results for the density 
of states at CNP are consistent in the limit of no Zee- 
man interactions and a perpendicular magnetic field (for 
6 = 0) with the experimental as well as theoretical results 
of [34-36]. 

Our results for both the Hall conductivity and the 
collisional conductivity are relevant to transport mea- 
surements performed at the CNP on epitaxial graphene 
grown on substrates such as SiC or BN where a band 
gap arises as a result of interaction with the substrate. 
In this regard, we have shown that one possible source 
of the opening of the gap in the density of states of 
the zeroth Landau level is the pseudo-Zeeman interac- 
tion which leads to the observed behavior of a plateau 
in the Hall conductivity and a dip in the collisional con- 



ductivity at CNP. In gaplcss graphene, such as graphene 
on Si02 substrate, opening of a gap in the density of 
states of the zeroth Landau level can also occur due to 
valley splitting of the Landau levels. In this case it can 
occur due to the inherent crystallographic symmetries of 
graphene [44]. The main difference in terms of realiza- 
tion in real physical systems is that in gapped graphene 
the effective Bohr magneton or the effective dipole mo- 
ment that couples with the extenal magnetic field is much 
larger which leads to larger valley splitting compared to 
gapless graphene for the same magnetic field. For an en- 
ergy gap of 0.28eV for graphene on SiC substrate, the 
effective dipole moment is about 30 times larger than 
the free electron spin magnetic moment. Therefore far 
smaller magnetic fields are required to observe the valley 
splitting in gapped graphene as against gaplcss graphene. 
V. SUMMARY 



In this work, we have investigated the coupling of an 
external magnetic field with the valley pseudo spin of 
Dirac fermions and its effects on electron transport in 
gapped graphene. Specifically, we have analyzed the 
splitting of the zeroth LL due to this pseudo-Zeeman 
interaction and its effects on the collisional and Hall 
conductivity at the CNP. To understand the role of the 
pseudo-Zeeman interaction we have obtained analytic ex- 
pressions and have plotted the results for the density of 
states at CNP in the presence of an external magnetic 
field whose tilt angle is varied. These results show that 
the pseudo-Zeeman interaction causes splitting of the ze- 
roth LL which vanishes when the magnetic field is aligned 
along the graphene plane. We find that the collisional 
conductivity at CNP shows activated behavior when the 
pseudo-Zeeman interaction is taken into account. Fur- 
thermore, we are able to show that as the temperature 
is increased for a fixed magnetic field, the closing of the 
gap in the zeroth LL occurs. 



VI. APPENDIX 



Equation (21) is written as 



_ 2{c lC2 ) 2 He 2 V F 



E 



(A.l) 



VM + (t z A z ) 2 /2H 2 lo 2 d - sVKI + (T' z A z )*/2hW D 
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For s, s' = +, + and +, — , the summation in the above equation for n = n + 1 is written as 



(/+*n — f-'n+l) 



(A.2) 



Note that n = n — 1 contribution vanishes as this corresponds to transition to filled states. Equation (A.2) can be 
simplified to yield 



E 



(^ + (^) 2 + V n + 1 + ^ )2 ) C/E»-/&h-i) 



™*i { + (■ V /?l + (^ 2 -V / " + 1 + (^ 2 ) - J 



(A.3) 



/ 



M> + (^) 2 + V / " + 1 + ^) 2 



For s, s = — , + and s,s = — , — , the summation on the right hand side of equation (A.l) is expressed as 

_ y^ (f-,n ~ 



(- V / " + (^7l) 2 -V / " + 1 + ^ 2 ) 
(/-% ~ 

Equation (A. 4) can be simplified to yield 

(- f + (^) 2 + V" + 1 + (^) 2 ) 2 - /Two 



(A.4) 



E 



+ (-V n+ (^ 2 -V n + 1 + (^ 2 ) ^-^i) 



(A.5) 



Considering the numerator of Eq.(A.3), we obtain 



r z A 



A 2 



(A.6) 



Z A 2 



2n+( z V ) 2 + ( z V ) 2 + 1 + 2Jn + ( ' * ) 2 Jn + l + ( 



'A, 

2 z \2 



X (f+,n ~ f+,n+l) 



x (/+% /-*n+l) 
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and then simplifying the numerator of Eq. 5, we get 



A, 



= -Jn + (^ + 1 n + l + (-^)2 (/::„- /T* +1 ) 



(A.7) 



+ - Jn + (-^=) 2 - Wn + 1 + (-^=) 2 (/:-„ - /:* +1 ) 



A 



*A 2 



:a 2 



2n+( z * y + ( g * Y + 1 - 2Jn + ( - > ) 2 W" + 1 + ( * % ) 



x (/-% — /+*„+! ) 



1 FujJdVZ 



x (/-*n — /-*n+l) 

The denominator of Eq.(A.3) and Eq.(A.5) is 



(^) 2 -y^+ 1 + (^5) a )' 



(r 2 - rl 2 - l) 2 . 



(A. 



One may notice that grouping terms such as +, + and +, — for s and s' that contain fl_" n leads to the cancellation 



of the following factor (2«/n + { ^^ ?\j n + 1 + (^f^Tf) 2 ) in Ec l- ( A - 6 )- Thc 

same holds for the — , — and 

terms in Eq. (A.7). Now using Eqs. (A. 3), (A. 5), (A. 6), (A.7) and (A.8) in Eq.(A.l), we arrive at the result 



(c lC2 ) 2 e 2 
ayx ~ nti X 



E 4 



n + 2 ( ti^D^) 2 + 2 ( fiJoV^) 2 + 2 
(/+% — /+%+! + — f-,n+l) 



/(r*-rf-l) 2 . 



Performing the summation over r z and t' z for +, + and +. — respectively, Eq. (A. 9) is simplified to yield 



(c lC2 ) 2 e 2 



E 4 

n 

E 4 



" + 2 ( fi j^V2 ^ + 2 ( fiJ^Va" ^ + 2 
(/+,n — + — f—,n+l) 

(/+,n — /+,n+l + /— ,n — f-,n+V 



> ■ 



After simplifying thc above equation we get 

\2 e 2 



( Cl c 2 ) 2 e 2 

^ " TTh X 



E 4 



-+(i±^) 2 + ^)x 



((2/+,n — /+,n+l + — /_, n+ i) — (/+,n + l + /-,n+l)) 



Similarly performing the summation over t 2 and t' z for — , + and — , — , Eq. (A. 9) is simplified to yield 



a ** = —^h~ x E 4 



"+^ + ) 2 +2-)x_ 
((2/+,n — /+,n+l + — f_ >n+ i) — {f+.n+i + /-,n + l)) 



Finally combining Eqs. (A. 11 and A. 12), we arrive at the final result for the Hall conductivity 



(A.9) 



(A.10) 



(All) 



(A.12) 



_ 2( Cl c 2 ) 2 e 2 

° vx ~ Vh x 



E 4 



x 

(/+,« — /+,»+! + f—,n — f—,n+l) 



with j = ±1. In the limit A z = 0, the above result reduces to that of [39] exactly. 



(A.13) 
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